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ABSTRACT 

The distribution of galaxy morphological types is a key test for models of galaxy formation 
and evolution, providing strong constraints on the relative contribution of different physical 
processes responsible for the growth of the spheroidal components. In this paper, we make 
use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupt¬ 
ing galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions 
usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., 
based on results from high-resolution hydrodynamical simulations, and we show that these 
new implementations reduce the efficiency of bulge formation through mergers. In addition, 
we compare our model results with a variety of observational measurements of the fraction of 
spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present 
uncertainties in the data represent an important limitation to our understanding of spheroid 
formation. Our results indicate that the main tension between theoretical models and observa¬ 
tions does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather 
from the distribution of galaxies of different morphological types, as a function of their stellar 
mass. 
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1 INTRODUCTION 

The morphological classification of galaxies constitutes one of the 
earliest attempts to explain the different observati onal propertie s 
of galaxies in terms of an evolutionary sequence ( lHubblelll92i) . 
The distribution of galaxies in different morphological classes has 


e.s. Dresslei 

1980l; 

Wilman et al. 20091) and stellar mass (see e.g. 

IVulcani et al 


Wilman & Erwin||2012|). Therefore modern as- 


trophysical research focuses on the complex interplay of physi¬ 
cal mechanisms acting on the baryonic components and respon¬ 
sible for the morphological transfor mation as a function o f stellar 
mass, cosmic epoch and environment terennan et alj201^ . Differ¬ 
ent approaches to the morphological classification of galaxy sam¬ 
ples have been proposed, moving from the visual classification of 
(relatively) small samples of local galaxies, to the use of automatic 
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pipelines applied to galaxy surveys. All these techniques have their 
strengths and shortcomings: typically the merit of the classification 
heavily depends on photometric data quality and/or on the algo¬ 
rithm used (e.g. profile or 2D fitting) and on the galactic compo¬ 
nents assumed. The aim is to distinguish between a centrally con¬ 
centrated spheroidal structure (also defined as a “bulge”) and a disc¬ 
like component (which may also host spiral arms). This oversimpli¬ 
fied picture neglects the heterogeneity of the bulge population and 
the role played by bar-like st ructures in galaxy evolution (see e.g. 
iKormendv & Kennicutll200^ . nevertheless it is useful to define the 
mass ratio between the bulge mass and the total mass of the galaxy 
(the so called bulge-to-total ratio, BjT) as a primary discriminant 
between bulge- and disc-dominated galaxies. 

Galaxies which consist purely of a bulge (ellipticals) are not 
easy to identify unambiguously in observations, even in the very 
local Universe. They have kinematic properties varying from pure 
pressure-supported systems to systems with partial rotational sup¬ 
port similar to le nticular galaxies, just without the defining o uter 
disk component jRrainovic et al.ll2008l : lEmsellem et alj|201 ih . In 
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absence of kinematics, an accurate multi-component decomposi¬ 
tion requires well resolved imaging at high signal to noise, for 
which a 2-component (bulge-rdisk) m odel can be fit (see e.g. 
ISimard et alj|201 lllMendel et al.ll20l4) . providing parameters tied 
to the morphological information. Visual classifications - confirm¬ 
ing the absence of disk light and asymmetric features such as bars 
and spiral arms - provides an alternative approach to identify el¬ 
lipticals. Nonetheless, ellipticals are always difficult to distinguish 
from face-on lenticular galaxies, and visual classification requires a 
huge effort for the characterization of large samples. I n both cases, 
with the average B/T increasing to high mass (see e.g. lBluck et al.l 
l2014h . it becomes more difficult to distinguish pure elliptical galax¬ 
ies from the increasing population of intermediate B/T (and still 
fairly concentrated) galaxies. Therefore, depending upon the exact 
cuts used to select elliptical or (more inclusively) early-type galax¬ 
ies, a different sample of elliptical galaxies may be selected. 

In a recent paper, IWilman & ErwinI ( l20I2h presented a lo¬ 
cal galaxy catalogue, based on a visual classification of galaxy 
morphology. This catalogue (SDSSRC3) includes galaxies with 
M* > 10^°'® Mq and it is based upon the New York Univer¬ 
sity V alue Added Galaxy Catalogue (NYU-VAGC. iBlanton et al.l 
bOOSh who mat ched the SPSS DR4 (Slo an Digital Sky Survey 
Data Release 4, lAdelman-McCarthv etal. II 20061) to the Third Ref¬ 
erenc e Catalogue of Bright G alaxies (RC3, de Vaucouleurs et ^ 
Il99lh . lwilman & ErwinI OoTl) used this sample to study the frac¬ 
tion of galaxies of different morphological types as a function of 
stellar and parent Dark Matter (DM) halo mass. For central galax¬ 
ies, the fraction of (visually classified) elliptical galaxies is found 
to increase with stellar mass and parent DM mass (only for central 
galaxies). They also found a relatively small fraction of genuine 
elliptical galaxies at high stellar masses (M* > lO^^'^M©), with 
a large contribution of SOs (which never have B/T > 0.7) and 
Spiral galaxies. One of the main differences between this work and 
previous ones lies in the definition of the galaxy sample (relatively 
nearby and bright - they appear in the RC3 catalogue) and in the use 
of visual classification to separate SOs from Ellipticals. Their def¬ 
inition is thus both qualitatively and quantitatively different from 
alternative approaches to define samples of Early Types galaxies 
bas ed on fitting light profile s and/or automatic decomposition (see 
e.g. lHvde & Bemardil2009l for the full SDSS volume). 

The work by IWilman & ErwinI l l2012t) has been extended in 
IWilman et alj (12013 ) including a detailed comparison to predictions 
from a suite of semi-analytic models (SAMs hereafter) of galaxy 
formation and evolution implementing different prescriptions for 
bulge formation. They showed that theoretical predictions are able 
to reproduce the observed increase of the Ellipticals fraction /e 
with stellar/halo mass. However, theoretical fs is too high with re¬ 
spect to the estimates from the SDSSRC3, by factors of a few (2.6 
to 4.2 for the models considered in that paper). Moreover, mod¬ 
els are also able to reproduce the fraction of passive galaxies and 
the fraction of star forming disc galaxies (the separation between 
star forming and passive galaxies being set at a specific star for¬ 
mation rate sSFR = 10“^^yr). Therefore, models overproduce 
Elli ptical galaxies at the expense of passive SOs and Spiral galax¬ 
ies. |^iman|et^i] l l2013h suggested a reduced efficiency of bulge 
formation in mergers, via efficient stripping of satellite galaxies 
leading to reduced baryonic merger ratios, as a viable solution to 
this problem. 

_ Contrasting results have been presented recently: IPorter et"^ 

(I 2 OI 4 I 1 compared pre dictions of improved variants of the 
ISomerville et alj dTOOSll mo del with the fract i on of Early Type 
galaxies computed using the iHvde & Bemai^ ( l2009h SDSS sam¬ 


ple, and conclude that their models tend to underpredict the fraction 
of Early type galaxies, when bulges are produced only in mergers. 
In view of these contradicting claims, it becomes important to un¬ 
derstand the origin of the differences between observational deter¬ 
minations of the Early Type/Elliptical fractions. 

The efficiency of binary galaxy mergers in destroying discs 
and forming bulges has been revised, following numerical sim- 
ulations taking into account the role of the gas content (see e.e . 

iRobertson et al.|[2006l : iGovemato et alj2009l : iHonkins et al.ll20()^ . 

Consensus has grown that the presence of a gas component may 
lead to the survival of the original disc even in 1:1 mergers. How¬ 
ever, these studies were based on idealized merger configurations 
and did not cover a statistic al sample of cosmol ogical initial con¬ 
ditions. As a matter of fact. iMoster et alj ( 1201 ih repeated similiar 
simulations including both cosmological accretion and accretion 
from a reservoir of hot gas associated with the pa rent dark mat¬ 
ter ha lo (two ingredients that were not considered in lHopkins et al.l 
l2009li and they did not find any clear dependence of burst efficiency 
on the progenitor’s cold gas fraction. 

A complementary approach has been used in t he framework 
of th e so-called Simulated Merger Tree Approach dMoster et al.l 
l2014h : this method is based on high-resolution hydrodynamical 
simulations of galaxy mergers, with initial conditions (Dark Mat¬ 
ter structures, orbital parameters, galaxy properties) extracted from 
(lower resolution) cosmological simulations coupled with SAMs. 
This technique thus combines the potential of high resolution to 
explore the effect of mergers on the mass distribution of the rem¬ 
nant galaxy, with cosmologically moti vated initial conditions. Us¬ 
ing this method, iKannan eTa D ( l2015h analyzed the channels of 
bulge growth in 19 simulations of binary galaxy mergers. In par¬ 
ticular, they followed the transfer of stellar and cold gas mass from 
the satellite to the central galaxy, as well as the mass exchanges 
between different components (disc, bulge, halo) of the central 
galaxy. They then try to quantify the mass flows by providing fitting 
formulae that parametrize the dependence of such mass exchanges 
on the merger mass ratio (both baryonic and in Dark Matter), and 
the orbital parameters. Their results show that the morphology of 
merger remnant depends mainly on the baryonic mass ratio be¬ 
tween the two merging galaxies but also on the gas fraction of the 
main galaxy, while the orbital parameters (i.e. the eccentricity of 
the orbit) play a minor role. 

The aim of this work is to study the implications of the new 
fitting formulae when implemented within models that can create 
large statistical realizations of the galaxy population, and in par¬ 
ticular the predicted fraction of early types as function of stellar 
and halo mass. This paper is organized as follows. In Section 
we introduce the semi-analytic implementations we consider in our 
analysis. We then compare the different distribution of morphologi¬ 
cal types to different sets of observational data in Sectionj^ Finally, 
we discuss our conclusions in Section|4] 


2 MODELS 

In this work we consider predictions from three different state- 
of-the-art SAMs, namely MORGANA dMonaco et alj l2007h . the 
IPe Lucia & Blaiz^ ( |2007| . DEB07 hereafter) version of the Mu- 
nich model and the Sant a Cruz model (SC-SAM . ISomerville et al.l 
l2008l : l^rter et alj|2014 1. In partic ular, we consider the MORGANA 
run defined in lLo Earo et al.l ( 120091 ) and calibrated on WMAP3 cos¬ 
mology (Ho = 72 km s“^ Mpc“^ , Clm = 0.24, ag = 0.8, 
n = 0.96, Da = 0.76); the DLB07 model applied to the Millen- 
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Figure 1. Bimodal distribution of predicted galaxy morphologies as a function of stellar mass. Contour levels mark number density levels (normalized to the 
maximum density) as i n the legend and the ho rizontal dotted line refers to B/T = 0.7. Filled coloured circles show the “classical” bulge to total ratios of 
observed galaxies in the lFisher & DrorvI j201 ih sample (i.e. we consider B/T > 0 only for those bulges classified as “classical” in the original sample). 


nium Simulation jSpringel et alj[200^ Hq = 73 km Mpc“^ 
, ~ 0.25, (Ts = 0.9, n = 1, = 0. 75) and the “fiducial” 

SC-SAM model as presented i nIPorter et alJ ll2014L inc luding also 
the modifications discus sed in Hirschmann et alj|2012l) applied to 


t alj 


20 111 Ho = 70.1 km s“ 


the Bolshoi simulation dKlypinet . 

Mpc~^ , nm = 0.279, 0-8 = 0.817, n = 0.96, Da = 0.721). 
We do not expect the differences in the cosmological parameters 
to affect model predict ions in a relevant w ay (see e .g Gu o et al.l 
l2013l;IWangetani2008l) . All models adopt a lChabrieJ ( l2003l) IMF. 
We refer the interested reader to the original papers for a com¬ 


plete overview of the modelling of the different physical processes 
in each code. Predictions from the same SAMs considered here 
have already been compared against a numbe r of observational 
constraints (see e.g iFontanot et alJl2009L [^1 ih . In the following, 
we give a brief description of the physical mechanisms relevant 
for bulge formation and morphological transformations. A detailed 
analysis of the timescales and channels of b ulge formation in MOR ¬ 
GANA and DLB07 h as b een presented in De Lucia et alj (1201 ll) . 
IFontanot et al.l | |201 ll) and lWilman et al.l ( l2013l) . 
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2.1 Reference and HOP09 runs 

All reference models assume that star formation happens primarily 
in disc systems. Spheroidal systems form when the stars and/or gas 
in the galactic discs lose angular momentum as a consequence of 
dissipative processes like galaxy mergers and/or disc instabilities, 
and are funnelled towards the centre of the galaxy. Broadly speak¬ 
ing, galaxy mergers are supposed to lead to the formation of an 
Elliptical galaxy or of a “classical” bulge (with kinematic and pho- 
tometric properties similar to those of Elliptical galaxies, see e.g. 
iDavies & Illingworthll983h . while secular processes are associated 
with the formation of “pseudo” bulges (c haracterized by disc-like 
expon ential profiles or kinematics; see e.g. lKormendv & Kennicuti 
|2004 . 

The modelling of disc instabilities represents one of the largest 
uncertainties in SAMs. Different assumptions are found in the lit- 
erature, from transferring just enough mass to restore stability (e.g. 
IPe Lucia & Blaiz^l2007h to the cat astrophic destructi on of the en¬ 
tire disc into a spheroidal remnant ( iBower et al.l l200 dl. Numerical 
simul ations of bar formation and evolution (see e.g. lPebattista et ahl 
l2006l) have not yet provided definitive evidence in favour of a given 
treatment or alternative prescriptions (although less dramatic mass 
transfers seem to be preferred Moster et al., in preparation); given 
these uncertainties, in the following we switch off disc instabili¬ 
ties in our reference models, that therefore correspond to the pure 
merger realizations in iwilman et al.l ( l2013h . It should be borne in 
mind that the contribution of disc instabilities is negligible at the 
highest stellar masses, but relevant for intermediate mass galax- 
ies (10^° < Mi,/M p) < 10^^) for all rn odels (as discussed in 
IPe Lucia et alj[20lll and IPorter et alJlioidh . thus playing a criti- 
cal role for the interp retation of observational data at these scales 
jFontanot et alj20Tll) . Moreover, in our models we also do not con¬ 
sider tidal stripping of stars in satellite galaxies. 

All 3 SAMs considered in this work implemented merger 
prescriptions based on a broad division between major and mi¬ 
nor mergers: in a major merger (i.e. mass ratio >0.3) the entire 
stellar and gaseous content of the two merging galaxies is given 
to the spheroidal remnant, which becomes an Elliptical galaxy 
(B jT = 1). In minor mergers, the stellar mass of the satellite 
is given to the bulge of the central/remnant, while the gaseous 
component is given either to the disc (DLB07, SC-SAM ) or to 
the b ulge (MORGANA ): i n both cases the merger triggers a star- 
burst. IWilman et alj (1201 3h showed that this scheme implies a sig¬ 
nificant increase of /b with stellar mass and parent halo mass for 
central galaxies. We refer to these model realizations as reference 
runs. In particular, for the SC-SAM reference run we use the same 
merger prescriptions as in I Somerville et al] l l2008l) . which includes 
the scattering of of a fixed (0.2) fraction of disc stars into the diffuse 
stellar halo (DSH) at each galaxy merger. 

_ A different formulation was proposed by iHopkins et al.l 

( l2009l) . based on resul ts from idealized hydr odynamical simula¬ 
tions of binary mergers teobertson et alj2006li . This “gas-fraction- 
dependent merger model” assumes that galaxy mergers trigger a 
burst of star formation, but the fraction of cold gas that participates 
in the starburst (i.e. efficiency of the conversion of cold gas into 
stars) depends on both the baryonic merger ratio of the two pro¬ 
genitors, and on the gas fraction of the disc of the primary galaxy. 
All stars formed in the burst are assigned to the spheroidal com¬ 
ponent of the remnant galaxy, as well as the whole stellar mass of 
the satellite galaxy. In addition, the coalescence of the two galax¬ 
ies might destroy a large fraction of the primary disc (depend¬ 
ing on the mass ratio of the merger), transferring its mass to the 


spheroidal component of the remnant. This is significantly differ¬ 
ent with respect to the reference approach, where the disc of the 
primary galaxy is either una ffected in minor merge rs or completely 
destroyed i n major merge rs. De Lucia et al.l j201lh studied the im¬ 
pact of the lHopldns et alj (l2009h prescriptions on bulge formation 
channels in the framework of the MORGANA and DLB07 models, 
finding small differences in terms of galaxy distribution into differ¬ 
ent morphological types, w ith respect to the reference runs. In con¬ 
trast, iHopkins et alj (l2009h have investigated the impact of these 
prescriptions in the SC-SAM model, finding a suppression of bulge 
formation in low-mass galaxies, and claiming that the gas depen¬ 
dence of merger-induced starbursts is a fundamental ingredient to 
reproduce the observed morphology-mass r elation. An improved 
version of the SC-SAM model inclu ding the Ho pkins et al.l l l2009li 
prescriptions has been presented in IPorter et~ ( l2014h . We refer 
to this set of realizations as the HOP09 runs. While running the 
reference and KAN15/Kz<l realizations for the SC-SAM model, 
we fixed all the oth er relevant parameters to the values assumed in 
IPorter et"^ ( l2014h . 

The frequency of mergers is clearly related to the timescales 
assumed for the orbital decay of satellites inside DM haloes. MOR¬ 
GANA and SC-SAM assign the merger time at the last time the 
galaxy is central (i.e. when the haloes merge), while in DLB07 
the residual merging time is estimated from the relative orbit of 
the two merging objects, at the ti me of subhalo disrup tion. MOR¬ 
GANA uses the fitting formula from iTaffoni et al.l l l2003h while both 
DLB07 and SC-SAM use variants of the dynamical friction for- 
mula. In particular, SC-SAM uses the formulation proposed by 
iBovlan-Kolchin et al.l (l2008h and DLB07 use the Chandrasekhar 
formula with a fudge factor 2 (which should bring merger times in 
close r agreement with those predicted by the lBovlan-Kolchin et al.l 
l2008l formula). The three approaches provide different pre dictions 
for the dynamics of satellite galaxies dPe Lucia et al .l201(]h . In par¬ 
ticular, the merging timescale for massive satellites (i.e. with large 
progenitor mass ratios) is roughly one order of magnitude shorter 
in MORGANA than in the other two models. 


2.2 Fitting formulae from iKannan et all d2015h . 

We briefly recall here the main results from iKannan et all (l2015l) . 
Using hydrodyna mical simulations base d on the parallel TreeSPH- 
code GADGET-2 dSoringel et alj|200^ . they show that in binary 
galaxy mergers, the mass exchanges between the satellite and the 
central galaxy and between the different components of the cen¬ 
tral galaxy primarly depend on the baryonic merger ratio (/ib) be¬ 
tween the two objects (computed at the infall time of the satellite 
galaxy, i.e. when it first became a satellite) and on the gas frac - 
tion (/gas) of the primary galaxy. Following iMoster et al.l d20l4) . 
initial conditions for the hydrodynamical simulations are set up at 
z ^ 1 using MORGANA to populate DM merger trees constructed 
from a cosmological volu me using the Lagrangian code PINOC- 
CHIO dMonaco et alj200j) . MORGANA provides predictions for the 
physical properties of merging galaxies, including stellar and cold 
gas mass, morphology and scale radii, while infall times are derived 
directly from the DM merger tree. It is wor th stressing that the sim¬ 
ulation suite used in iKannan et al.l d2015h analysis includes only 
19 merger configurations, extracted from the cosmological volume 
and bound to have DM halo masses at 2 = 0 similar to the Milky 
Way host halo. Initial masses for primary galaxies lie in the range 
9.08 < log(M*) < 10.77, with merger ratios 0.01 < /ru < 1 and 
gas fraction 0.1 < /gas < 0.6. This is a rather small sample and 
results in a significant scatter around the proposed analytic fitting 
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Figure 2. Elliptical galaxy fraction (B/T > 0.7) as a function of stellar mass {upper panels) and parent halo mass {lower panels). Upper and lower rows 
refer to central and satellite g alaxies respectively. In each panel, datapoints with eiTorbars con'espond to visually classified total derived from the SDSSRC3 
sample iWilman et alj|201~3l . Solid, dot-dashed, dashed and dotted lines refer to the SAMs predictions for the reference, HOP09, KAN15 and Kz<l runs, 
respectively. Only bins with more than 10 model galaxies have been considered. 


formulae, thus limiting the precision of the resulting fitting formu¬ 
lae. Moreover, since we will apply these prescriptions on a full cos¬ 
mological volume, we must extrapolate beyond the DM halo mass 
scale and redshift range they were calibrated on. Nonetheless, given 
the present uncertainties in the modeling of bulge formation during 
galaxy mergers, it is interesting to explore the implications of these 
trends in terms of the distribution of model galaxies in the differ¬ 
ent morphological types, in order to assess if these results go into 
the direction of relieving or exacerbating the tension between data 


and model predictions. In order to test the strength of our results, 
we define additional model realizations imposing a (gaussian dis¬ 
tributed) scatter around the adopted fitting formulae corresponding 
to variance ranging from 0.2 to 1. The analysis of these runs show 
that none of the conclusions we discuss in this paper change when 
scatter in the fitting formulae is taken into consideration. 

In detai l , the main mass exchanges parametrized in the 
iKannan et alj (1201 Sh simulations include the fraction of total (stel- 
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Figure 3. Galaxy stellar mass functions at 2 : = 0. Line types and colours refer to SAM predictions as in fig. |2] Light grey points refer to the stellar mass 
function compilation from iFontanot et al] j2009l see references herein). 

lar plus cold gas) satellite mass ending up in the bulge of the rem¬ 
nant galaxy 

/S2B = Mb, (1) 

where 1 —/s 2 b represents the satellite mass deposited into the DSH 
as a result of tidal stripping acting on the secondary galaxy, while it 
orbits within the parent DM halo potential well. They also consider 
the fraction of the central disc gas mass transferred to the central 
bulge 

/g 2B = (1 — /gas)Mb, (2) 

which is strictly defined as < 1. The cold gas transferred from the 

primary disc to the bulge of the remnant, and the cold gas from the 
satellite are readily consumed in a burst of star form ation, follow¬ 
ing the original definition in each model. Finally, the lKannan et al.l 
(l2015h fitting formulae describe the mass loss from the stellar disc 
of the central galaxy, as the result of tidal effects: i.e. the fraction 
of primary disc stars transferred to the bulge of the remnant 

/d2b = 0.37Mb, (3) 

and the fraction of central disc stars dispersed into the DSH 

/d2h = 0.22Mb. (4) 

It is worth stressing that, for massive galaxies, results are more sen¬ 
sitive to the fitting formulae describing the effect of mergers on 
existing stars, as they typically have very low gas fractions. 

In the following, we focus on two different sets of models: we 
first co nsider SAMs realizations (KAN15) where the lKannan et al.l 
(l2015h fitting formulae are implemented at all redshifts. In a second 
set of modified models (Kz<l), we allow these modifications only 
for z < 1, keeping the reference merger prescriptions at earlier 
times. This second choice is motivated by the fact that the orig¬ 
inal simulations considered in iKannan et al.l ( 1201 5h are set up at 
z = 1 (using the information provided by the SAM) and then run 


to z = 0. Since we expect the properties of higher-redshift discs to 
be different than those of low-z systems , Kz<l represents a con¬ 
servative lower limit on the effect of the lKannan et al.l j2015l) pre¬ 
scriptions on the distribution of morphological types. 


3 RESULTS & DISCUSSION 

In fig. [U we compare fhe bidimensional M* vs B/T distribu¬ 
tions in our model galaxy catalogues and in the local (closer 
than 11 Mpc) v o lume- limited sample of 321 nearby galaxies from 
iFisher & DrorvI (1201 ih . based on a bulge-disc decomposition and 
pseudo-bulge diagnostics in the Spilzer 3.6Mm band. As disk insta¬ 
bility has been switched off in the models, we only consider “clas- 
sic al” B/T ratios for the data sample (i.e. we consider B jT > 0 in 
the lFisher & DroijgOlT sample only for those bulges classified as 
classical, see also lWilman et alj|2013l) . A consistent picture can be 
drawn from the predictions of the different reference models, de¬ 
spite different merger prescriptions and merger timescales. In the 
KAN 15 runs (i.e. when the new prescriptions are applied at all red- 
shifts) the abundance of B/T >0.7 galaxies is drastically reduced 
at all mass scales. However, this is obtained at the expense of almost 
completely devoiding the B/T > 0.9 region. This effect goes into 
the same direction as for the HOP09 runs, but the ov erall reduction 
of bulg e dominated galaxies is greater when using the lKannan et al.l 
( l2015h prescriptions in these models: the KAN15 column of Fig.[T] 
shows that most model galaxies have B/T < 0.5. As expected, the 
Kz<l model provides intermediate results between the reference 
and KAN15 runs: the number of B/T > 0.7 galaxies is reduced, 
but not as much as in the KAN15 runs. In particular the effect is 
mass dependent, as more massive galaxies are scattered to a wider 
distribution of B/T with respect to the reference runs (where they 
are mostly B/T = 1 objects). 

More insight can be obtained by comparing the predictions of 
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our models with the observed f g as a function of ste llar and par¬ 
ent DM halo mass, as defined in iwilman et alj l l2013tl . As in their 
work, in Fig. [21 we define as Ellipticals all model galaxies with 
B jT > 0.7, and split our sample by hierarchy considering sep¬ 
arately central and satellite galaxies in the upper and lower rows 
(only bins containing more that 10 objects are considered). Our 
choice for the threshold at B/T ^ 0.7 is motivated by ob serva- 
tional studies dWeinzirl et alJlioO^ : iLaurikainen et alj|201oh ) and 
by consistency with previous studies. We want to stress that our 
results do not change qualitatively if we assume a different B/T 


threshold for the definition of Elli ptical galaxies (in the ra nge 0.5- 
0.9). We remind the reader that the lWilman & ErwinI ( 1201 2h sample 
contains 854 galaxies with M* > IO^^’^Mq, whose morphologies 
are visually classified. It is difficu lt to directly compare this sam¬ 
ple with the lFisher & Dror^ 1 I2OIII) catalogue, as the latter contains 
only 11 galaxies with M* > 10^° ®Mq, while the former does not 
allow for a finer classification between pseudo and classical bulges. 
It is thus worth stressing that in the following plots fs formally 
represents an upper limit for the fraction of Elliptical galaxies, as it 
is not possible to discriminate the relative contribution of disk in¬ 
stabilities and merger proces ses in the formatio n of a spheroid, but 
only its main shape (but see lErwin et'^l2015l for recent progress 
on these issues). Nonetheless, a more detailed classification should 
have a limited impact on our conclusions, since we do not expect 
the formation of B/T >0.7 objects to be driven by disk insta¬ 
bilities alone, although they might contribute to increase the mass 
of a merger-induced spheroid. Similar arg uments hold for the other 
two main samples c onsidered in this work 1 H yde & Bemardil20091 : 
[Mendel et al.ll201^ . As discussed in iwilman et al.l2013l) . the ref- 
erence models (solid lines) overpredict /_b over the whole mass 
range. The discrepancy is larger for more massive systems, where 
model fractions almost reach 100%. The discrepancy is particu¬ 
larly severe for central galaxies that show, in all models, a strong 
dependence of fs on stellar and parent halo mass. It is interesting 
to notice that the refer ence run of SC-SAM (that was not included 
in iwilman et al.l2013l) predicts lower elliptical fractions at the high 
mass end, with resp ect to the other two m odels, thus providing the 
closest match to the lWilman et al.l l l2013h sample among reference 
runs. Moreover, the predicted fractions of B/T > 0.7 galaxies in 
the HOP09 runs (dot-dashed lines) deviate considerably from the 
reference runs only for the SC-SAM model, while they are very 
close for the other two models (as already noticed in IWilman et al.l 
l2013h . 


The picture changes when the KAN 15 runs are considered; 
the predicted fs are reduced at all mass scales, both as a func¬ 
tion of stellar and parent halo mass. It is worth stressing that the 
final masses of model galaxies change between the reference and 
the KAN15/Kz<l runs (dashed and dotted lines), due to the dif¬ 
ferent assembly history. The effect of these changes is imprinted 
on the z = 0 galaxy stellar mass function (Fig. [3}. For both MOR¬ 
GANA and DLB07 there is a clear decrease in stellar mass at the 
high-mass end due to the deposit of stellar mass into the DSH as a 
consequence of the assumed modeling for fs2B and fD2H (and be¬ 
cause more massive galaxies experience more merger events). This 
effect is more severe for DLB07, as the number density of massive 
galaxies is already the smallest among the models considered and 
the original prescription of this model did not include stellar strip¬ 
ping of satellites (and we did not retune the parameters to account 
for this effect). In this model, both the KAN15 and Kz<l runs do 
not predict any galaxy with stellar mass M* > IO^^'^Mq, and the 
highest mass bins available are affected by small number statistics. 
In addition, this model shows the smallest reduction of /b, with 


the M* > Mq bins still dominated by B/T > 0.7 galaxies. 
We checked that, in DLB07, this behaviour is mainly due to the re¬ 
duced mass of infalling satellites (i.e. /s 2 s). We also noticed that 
the HOP09 runs of MORGANA and DBL07 predict galaxy stellar 
mass functions almost indistinguishable from those of the reference 
runs. For the SC-SAM instead, only small differences between the 
predicted mass function from the four runs can be seen. This is 
due to the fact that the SC-SAM already inclu des the scattering 
of disc star s into the D S H at e ach merger in the [Somerville et aP 
(l2008i) and Porter et alj l l2014tl versions. Therefore, the impact of 
the Kannan et alj 1 20151) formulae is reduced, with respect to the 
other two models. 

Despite for MORGANA and SC-SAM the predicted frac- 
tions for the KAN 15 r uns are now roughly consistent with the 
IWilman & ErwinI l l2012tl observational constraints, this is obtained 
at the expenses of losing all B/T > 0.9 Ellipticals. The Kz<l 
runs provide intermediate results with respect to the other two sets 
of runs. In most cases, h owever, the predicted /e still overpredicts 
the observed fs in the Iwilman & ErwinI ( 1201 2l) sample, at most 
mass scales. 

We then compare our model predictions with additional ob¬ 
servational constraints in Fig.[4| We first consider different B/T 
ranges for model galaxies: we split the sample into B/T > 0.7 
(upper row) and B/T < 0.7 objects (mid row), and we also con¬ 
sider a B/T <0.1 selection (lower row), which corresponds to the 
fraction of “bulgeless” galaxies, i.e. those objects with the smallest 
contribution from a spheroidal component. As for the observational 
data, we consider four different samples and w e plot their results 
in eac h panel as follows. We first consider the IWilman & ErwinI 
( l2012h sample, and we plot as solid circles their total /e in the 
upper row. In the mid row, we plot the sum of the (visually clas¬ 
sified) SO and Spiral fractions. We include as open squares in the 
upper row, the total fraction of Ear ly Type galaxies in the SDSS 
selected bv iHvde & Bernard ( l2009ll on the basis of the SDSS pa¬ 
rameter fracDev=l (the fraction of source light well-fitted by a de 
Vaucouleur’s profile) and r-band axis ratio b/a > 0.6. 

We also consider the catalogue of bulge, disk and total stel¬ 
lar mass est imates for ga l axies in the Legacy area of the SDSS 
presented in iMendel et al.l ( l2014ll . Those mass estimates are based 
on updated SDSS photom etry and bulge-(-disc decomposition, us- 
ing dual Sersic p rofiles dSimard et alj 1201 ih . In their analysis, 
iBluck et^ ( I 2 OI 4 I 1 consider a subsample of the full catalogue de¬ 
fined by the cuts 0.02 < 2 < 0.2, 8 < log(Mstar) < 12, 
dDB < 1 and Pps < 0.5. dDB represents the scatter in the rela¬ 
tion between the estimate of the total stellar mass and the sum of the 
estimates of bulge and disc separately; the applied cut removes ob¬ 
jects with problematic decomposition. Pps is the probability that a 
source is best-fit by a single Sersic profile; the applied cut removes 
sources with possibly spurious compon ents. It is wor t h stre ssing 
that the bulge-(-disk modeling adopted in iMendel et al.l ( |2014|) . i.e. 
the combination of a de Vaucouleur’s bulge and an exponential 
disk profiles, is unable to accurately recover the total flux of galax¬ 
ies with Sersic index larger than 5. For these sources, the model 
fails in reproducing the bright inner and extended outer profiles, at 
the same time. Nonetheless, we include also these objects in our 
analysis, as they might represent some extreme configurations of 
bulge/disc-dominated sources. In order to assess for the uncertain¬ 
ties i n their bulge-(-disc de composition, we define two subsamples 
of the lMendel et al.l ( l2014l) catalogues as f ollows. We first con sider 
all galaxies satisfying the first 3 cuts as in iBluck et alj ( l2014h . We 
use this catalogue to define a subsample (A) imposing B/T = 0 
to all objects with Pps ^ 0.5, in order to obtain a lower limit 
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MORGANA DLB07 SC-SAM 
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Figure 4. Distribution of galaxies among morphological types as a function of stellar mass. Filled circles refer t o the t otal fraction of Ellipticals in iwilman et alJ 
(20131), while open squa res to the total fraction of Early Type galaxies from the sample_of|Hj;de_&Bernard (20091) . Shaded areas represent the results from 
the iMendel et alj (2oi4h decomposition, while stars and hatched histograms refer to lEisher ifePronTIlonl) sample. See main text for more details on the 
observational samples considered. Line types and colours refer to SAM predictions as in fig. Right-hand labels indicate the model galaxy BjT range 
considered in each row. 


on the fraction of bulge-dominated galaxies in the iMendel et al.l 
) l2014h sample. Analogously, we define a subsample (B) assum¬ 
ing B/T = 1 for sources with Pps ^ 0.5, which maximizes the 
fraction of bulge-dominated galaxies. The shaded area in Fig. |4] 
represents the confidence region between the resulting fractions of 
B/T > 0.7 and B/T < 0.7 objects in the two subsamples. Frac¬ 
tions referring to subsample A correspond to the lower and upper 


envelope of the shaded region in the upper and mid panels respec¬ 
tively (viceversa for subsample B). 

Finally, we include results from iFisher & Dro^ (201 ih . This 
local sample allows a finer classification of spheroids into “classi¬ 
cal” and “pseudo”-bulges, thus providin g a finer comparison with 
our model predictions with respect to the lMendel et alJ (20141) sam¬ 
ple, but it contains only 11 galaxies in the stellar mass range of in¬ 
terest (i.e. Mi, > lO^^’^M©). As in all our runs we switched off 
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disc instabilities, our models formally account for the formation of 
classical bulges only, i.e. we assume that the entire “pseudo”-bulge 
population is the result of secular processes. We first plot the result¬ 
ing morphological fractions as stars in the upper and middle panels 
of Fig.|4] These fractions do not depend on the bulge classification, 
since there are no galaxies hosting pseudo -bulges with j3/T > 0.7 , 
and in reasonable agree ment with both IWilman & ErwinI l l2012tl 
and iMendel et al. ^20141). T he hatched histograms in the lower row 
represents (as in lpontanot e t alJl201 ih the fract i on of galaxies with 
classical B/T < 0.1 in the Fisher & DrorvI bOllh sample. We 
stress that the sample contains a significant number of galaxies 
hosting pseudo-bulges with B/T > 0.3, which may hide a (sub¬ 
dominant) merger-driven classical bulge. This effect would tend to 
reduce the estimated fraction of bulgeless galaxies and the hatched 
histograms then represent an upper limit to the bulgeless fraction in 
the sample. 

This comparison is meant to account for the uncertainty in 
the definition of galaxy morphology. In general, there are still 
large differences in the observational estimate s, which explain the 
app arently contrad i ctory conclusions given in IPorter et al.l l l2014ll 
and IWilman et alj ( l2013h . The tension between different sam¬ 
ples is likely due to the different techniques adopted to estimate 
the relative importance of bulges and discs. The biggest discrep- 
an cy is seen in the upper r ow between the Early Type catalogue 
of Hj^de^jenwdi|j 20091) and the Elliptical fraction from the 
IWilman & ErwinI l l2012h sample: the lMendel et al.l ( l2014h decom¬ 
position shows intermediate results, but is in better agreement with 
the latter sample at the high mass end. In the mid row, we show how 
the reduced bulge formation efficiency in galaxy mergers impacts 
the statistics of disc-dominated galaxies. As in Fig.[2 we see that 
the KAN15 run of the MORGANA model and the SC - SAM runs pro¬ 
vide the best agreement with the lWilman & Er\t^ l l2012tl sample, 
also when B/T < 0.7 model galaxies are compared with sources 
with a relevant disc component (either SOs or spiral galaxies). How¬ 
ever, none of the models and runs c onsidered in this wor k is able to 
reproduce the peculiar trends of the lMendel et al.l (12014h sample. 

Finally, in the lower row, we focus on the number density 
of “bulgeless” galaxies. The formation of these objects has been 
claimed to represent a pot ential challenge for th e standard cos¬ 
mological model (see e.g. Kormen^ et_alJ l201Ct and references 
herein). We showed in iFontanot et alj ( 1201 ll) that these objects do 
not constitute a rare population in SAMs, as there is a sufficient 
number of DM halos with merge r-quiet assembly h istories to ac¬ 
count for local statistics (see also lPorter et alj|2014l) . With respect 
to the reference runs, both the KAN 15 and Kz<l realizations pre¬ 
dict an increase of the fraction of B/T < 0.1 galaxies at all 
mass scales, which improves the agreement with the observations. 
The increase is relatively mild in DLB07 and SC-SAM and more 
marked in MORGANA . This is due to the shorter timescales for 
merging assumed by the latter model, which induce more galaxy 
mergers. 

A major role in the lKannan et alj ( l2015ll formulation is played 
by the stars getting stripped (from the satellite galaxies) or scattered 
(from the primary discs), and dispersed into the parent DM halo. In 
Fig. a we show the stellar mass in the DSH as a function of halo 
mass, normalized to the stellar mass of the central galaxy plus the 
DSH itself. The upper panel refers to KAN 15 runs, while the lower 
panel to Kz< 1 runs. In all cases, there is a clear trend of a larger 
contribution of DSH at increasing halo masses. A direct compari¬ 
son of these DSH estimates with observational constraints is com¬ 
plicated by the intrinsic difficulties in disentangling the contribu¬ 
tion to the total cluster luminosities from the central bright galaxy 



Figure 5. Fraction of stellar' mass scattered into the DSH component (rel¬ 
ative to the mass of the central galaxy plus the DSH) as a function of the 
mass of the parent DM halo. 


and from the DSH (i.e. the intra clus ter light), due to the different 
techniques and assumptions (see e.g. lZibettill2008l . for a review on 
these issues). At face value, our pr edictions agree well with con¬ 
straints from lOonzalez et alj (120(151 ^ 0.6 at ATdm ~ 10^ ^AT(7)) , 
but are larger than the recent estimates from iBudzvnski et al.l(l2014L 
~ 0.2 — 0.4). Moreover, the differences between the KAN15 and 
Kz<l runs are small, showing that in these models, most of the 
DSH is formed at z < 1. This confirms that most of the differ¬ 
ences at the high mass end between the reference models and the 
new runs are due to the effect of stellar stripping/dispersion, which 
reduce the mass of both progenitors, rather than to the influence 
of gas fraction (as noted above, most model massive galaxies have 
typically low gas fractions at z < 1). 


4 CONCLUSIONS 

In this paper we study the impact of a different modelling of bary- 
onic (stellar and cold gas) mass transfers associated with galaxy 
mergers on the distribution of morphological types as a function 
of hierarchy, stellar mass and parent DM mass. We consider three 
independently developed SAMs, namely the SC-SAM model, the 
DLB07 version of the Munich model and MORGANA . When com¬ 
pared to a local sample of visually classified galaxies, the two lat¬ 
ter models have already been shown to overpredict the fraction of 
Elliptical galaxies at high masses, while reprodu cing the observed 
fraction of passive galaxies. Iwilman et al.l ( l2013ll ascribed this dis¬ 
crepancy to a too high efficiency of morphological transforma¬ 
tion (from disc- to bulge-dominated system) during galaxy mergers 
(which they suggest can be reduced by assuming that a fraction 
of the satellite galaxy is deposi ted into the DSH b efore the merger 
takes place). On the other hand. lPorter et al.l ( l2014b showed that the 
SC-SAM model underpredicts the fraction of Early Type galaxies 
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in the lHvde & BemM^ ( l2009h SDSS sample. In this work, we ex¬ 
tend the analysis by com paring model predictions to four differ¬ 
ent observational sample s iHvde & Bernardil2009l.lu sher & DrorvI 
I2OI lHwilman&Erwinl [2012!: iMendelet all l2014h and by con 
sidering new prescriptions fo r the mass exchange s during galaxy 
merges, recently proposed bv iKannan et al.l ( l2015h . These authors 
analyzed high-resolution hydrodynamical simulations of 2 < 1 bi¬ 
nary mergers, with initial condition extracted from a combination 
of cosmological simulations and SAMs. The simulations were used 
to parametrize the strength of baryonic mass transfer between the 
different components (disc, bulge, halo) as a function of an effec¬ 
tive merger ratio (taking into account both the intrinsic merger ratio 
and the details of the relative orbit). 

Our results confirm that the different modelling of mass trans¬ 
fer during mergers leads to significant modification of the distri¬ 
bution of B/T ratios as a function of both stellar and parent halo 
mass. In particular, the fraction of B/T > 0.7 ga laxies, /g, is 
reduce d at all mass scales. However, runs using the lKannan et al.l 
(l2015h fitting formulae at all redshifts result in virtually no galaxy 
with B/T — 1. Moreover, in these model realizations most of the 
bulge mass is assembled at 2 > 1. More conservative runs, en¬ 
abling the new fitting formulae only at 2 < 1 (i.e. the redshift 
range probed by the simulations) provide a better match to the B/T 
distribution, at the expense of a smal ler reduction of /g. T hese re¬ 
sults suggest the need to extend the iKannan et alj ( 1201 Sh work at 
higher redshifts. In fact, 2 > 1 discs are lik ely very different struc¬ 
tures with respect to local Spirals (see e.g. iFdrster Schreiber et al.l 
l2009ll . therefore we expect them to show a different response to the 
tidal fields triggered by close encounters, with respect to their local 
counterparts. 

The detailed effects of the new prescriptions on the distribu¬ 
tion of galaxy morphologies are strongly SAM dependent. MOR¬ 
GANA shows a significant reduction of fe which moves the m odel 
predictions in better agreement with the lWilman et alj ( 1201 3h data. 
In DLB07, the most relevant effect is that of decreasing the mass of 
the most mas s ive ga laxies, due to the stripping effects included in 
iKannan et al.l(l2015l) formulae. At intermediate-to-low mass scales, 
model predictions are similar to the SDSSRC3 fractions, while 
at higher masses they are closer to the early type fraction in 
iHyde & Bemardil ( l2009ll than to the fs estimate from SDSSRC3. 
Finally, the effect of the new prescriptions is limited in SC-SAM 
due to the fact that the reference run already predicts a reduced fs 
with respect to the other two SAMs we co nsider: this model still un - 
depredicts the Farly Type fractions as in iHvde & Bernard ( l2009ll . 
but its predictions are similar to the fs values from SDSSRC3. 

Overall, the new formulation predicts a decrease of the frac¬ 
tion of bulge-dominated (B/T > 0.7) galaxies and a correspond¬ 
ing increase in the fraction of disc-dominated galaxies (including 
bulgeless galaxies). The comparison between model predictions 
and data is complicated by the different observational estimates 
for the fraction of bulge-dominated objects. As an additional test, 
in this work we also cons i der th e morphological fractions com¬ 
ing from the iMendel et al.l ( |2014|) sample. In general, this sample 
provides m orphological fractions w hich are intermediate with re - 
spect to the lWilman & Erwiiil ( 1201 2h and iHvde & Bernard j2009tl 
datasets, but at the high mass end, bulge-(-disk decomposition in 
this sample is consistent with the results based on visual inspec¬ 
tion of SDSSRC3 galaxies. It is worth stressing that the spread be¬ 
tween the different observational constraints for ]e is significant. 
Much of this tension is directly linked to the different techniques 
used to estimate galaxy morphology, and ranging from purely auto¬ 
matic classification to human visual inspection. Additional sources 


of uncertainty come from the presence of fine structure features like 
spiral arms and bars, whose detection strongly depend on the data 
quality the classification is performed on (independently from the 
overall strategy). Much work is therefore needed, on the observa¬ 
tional side, in order to assess the consistency of the different clas¬ 
sifications and/or to understand the systematics involved in each 
choice. The end product of this process will provide even stronger 
constraints for modelers trying to quantify the amount of morpho¬ 
logical transformation connected to the different channels of bulge 

growth. _ 

We stress that the IKannan et al. I |2 oT^ fitting formulae are 
based on a small sample of high-resolution hydrodynamical simu¬ 
lations (19 binary isolated mergers, covering a fair range of merger 
ratios), and their results show a large scatter around the fitted rela¬ 
tions. More high-resolution simulations are thus required to better 
quantify both the mean relations and their scatter, which could also 
play a role in defining the observed distribution of morphological 
types. Moreover, galaxy mergers in SAMs are assumed to be bi¬ 
nary, i.e. involving only two galaxies. This is an idealized assump¬ 
tion useful for most physical configurations, but not representative 
when more than one satellite interacts with the central galaxy at the 
same time (see e.g. lVillalobos et aDl2014l) . More work in this di¬ 
rection is needed to characterize mass transfers in multiple merger 
configurations, as a function of the relative masses and orbits. 
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